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Protein allostery is essential for vital activities. Allosteric regulation of human hemoglobin (HbA) with two 
quaternary states T and R has been a paradigm of allosteric structural regulation of proteins. It is widely 
accepted that oxygen molecules (O2) act as a "site-specific" homotropic effector, or the successive O2 
binding to the heme brings about the quaternary regulation. However, here we show that the site-specific 
allosteric effect is not necessarily only a unique mechanism of O2 allostery. Our simulation results revealed 
that the solution environment of high O2 partial pressure enhances the quaternary change from T to R 
without binding to the heme, suggesting an additional "non-site-specific" allosteric effect of O2. The latter 
effect should play a complementary role in the quaternary change by affecting the intersubunit contacts. 
This analysis must become a milestone in comprehensive understanding of the allosteric regulation of HbA 
from the molecular point of view. 

Complex processes of signaling in the living cell are comprised of many allosteric interactions between 
proteins \ Thus controlling the allostery is emerging as a promising approach for pharmaceutical drug 
discovery and allosteric ligands for G protein-coupled receptors have been approved as medicines^. 
Allosteric regulation is assumed to be triggered by site- specific ligand binding to allosteric sites, which are distinct 
from the orthosteric ligand binding site. As a paradigm of such allosteric regulation, allostery of human hemo- 
globin (HbA), which is an oxygen transport protein composed of two ot and two P subunits, has long been 
thoroughly studied^ Under the circumstances, it is widely accepted that O2 acts as a "site-specific" homotropic 
allosteric effector, or the successive O2 binding to the four subunit hemes triggers quaternary change from the low 
O2 affinity T state to the high affinity R one. The MWC modeP is one of the most successful models to explain the 
sigmoidal oxygen binding curve on the basis of the homotropic allosteric mechanism. From the molecular point 
of view, the homotropic allostery was additionally explained by including the interaction energies of intersubunit 
contacts which stabilize the T state structure^; when O2 binds to the heme, the tertiary structural change of the 
HbA subunit is triggered by the heme deformation from the domed structure to the planar one and, as a result, 
those contacts are destabilized to drive the T to R transition^ In fact, in HbA subunit analogous protein, 
myoglobin (Mb), it has been elucidated that the heme deformation by the ligand (O2 or CO) binding or 
dissociation triggers structural relaxation in the whole protein^" However, our current molecular dynamics 
(MD) simulation revealed that, even without site-specific O2 binding to the heme, the solution environment of 
high O2 partial pressure enhances the quaternary structural change from T to R, suggesting the "non- site- 
specific" allosteric effect of O2. 

Using the statistical ensemble MD approach^ we executed 128 independent 8 nanosecond (ns) MD 
simulations starting from the T state deoxy-HbA crystal structure^^ at the physiological condition (310 K and 

1 atm) in both aqueous solvent (MD"^^* with ~ 12000 H2O molecules) and oxygen-rich aqueous solvent with [O2] 
= 0.55 mol/L (MD°2 with 120 O2 and -12000 H2O molecules), which is about 500 times higher than that in the 

02 saturation concentration at ambient conditions, ~1 mmol/L. This high partial pressure was applied to 
amplify the non-site-specific effect of O2. We artificially ignored the O2 binding to the heme and thus all the 
hemes were retained in the domed structure for the purpose to exclude the heme O2 binding effect. 

Results 

Global evidence of quaternary change: transition in intersubunit dihedral angle and RMSD. It was reported 
that spontaneous T to R quaternary transition reproducibly occurred in the submicrosecond MD simulations of 
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Figure 1 | T state crystal (initial MD) structure of HbA (PDB ID: 2DN2). 

The dihedral angle x defined by the four centers of mass, i.e., those of oci 
subunit, ociPi dimer, a2|32 dimer, and 0^2 subunit, is shown as ball and stick 
drawing. The blue color drawing is the COMs of the initial T state structure 
with X = —99.5 degrees. The red one is those from a typical T to R 
quaternary change trajectory at 4.584 ns with x = —87.8 degrees. 

HbA starting from the T state structure^^'^^. Our simulations also 
exhibited the T to R quaternary changes. Because the quaternary 
change is accompanied with ajpi (a2p2) dimer rotation^'^^ the 
temporal change of the intersubunit dihedral angle defined by 
the four centers of mass (COMs) of subunits (see Fig. 1), is an 
index of the quaternary change; —99.5 and —89.0 degrees for T 
and R crystal structures (PDB ID: 2DN2 and 2DN1''), respectively. 
Figs. 2 gives one of the typical trajectories of MD°^ in which the T to R 
quaternary change occurred. Until 3 ns, the HbA structure kept its 
initial T structure, and later, x changed during t = 3.0 to 4.5 ns and 



reached —87 degrees, strongly correlating with the decrease of 
RMSD^ to 1.4 A (Fig. 2a). 

Local evidence of quaternary change: switch region rearrangement. 

Furthermore, local structural changes in the aiP2 (0^2 Pi) switch region 
were also consistent with the T to R change: P2His97 moves from the 
groove between aiPro44 and aiThr41 to that between aiThr38 and 
otiThr41^'^^. This rearrangement also occurred in this trajectory 
during t = 3 to 4.25 ns (Fig. 2b). The aiThr38-P2His97 distance 
decreased from 11 to 5 A and that of aiPro44-p2His97 conversely 
increased from 5 to 11 A (Fig. 2c). Other trajectories of MD^^ and 
jy[j-)WAT ^hich the quaternary change occurred also revealed the 
similar behaviors of the RMSD^, and switch region structural 
changes (Supplementary Figs. SI and S2). 

Other regions also revealed T to R quaternay change characteris- 
tics. Breakages of intersubunit salt bridges aiLys40-p2Hisl46 and 
aiTyr42-P2Asp99 occur by the T to R quaternary change^^. These 
breakages were also observed in the quaternary change trajectory 
(Supplementary Fig. S3). Meanwhile, in the otiP2 (ot2Pi) hinge region 
(also denoted as flexible joint region), the intersubunit contacts were 
retained and there were no significant structural rearrangements 
during the quaternary change (Supplementary Fig. S4) as crystal- 
lographically observed^^. 

Non-site-specific allostery emerged in the intersubunit dihedral 
angle distribution. To investigate the effect of high O2 partial 
pressure, statistical comparison between the 128 MD trajectories of 
MD°^ and MD"^""* was carried out. When the time average of x over 
7-8 ns is larger than —95 degrees, the trajectory is defined as a 
quaternary change trajectory. On this definition, the number of 
quaternary change trajectories were 10 and 4 for MD°^ and MD'^''^ 
respectively, suggesting that the high O2 partial pressure enhances 
the T to R quaternary change. 
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Figure 2 | T to R quaternary structural change. A typical quaternary change trajectory of MD*^^. (a) Changes of the RMSD to the T (RMSD^) and R 
(RMSD^) crystal structures and the intersubunit dihedral angle x- (b) Structural rearrangement at the switch region of the ai|32 interface. Important 
residues are drawn at 3 ns (red) and 4.25 ns (blue) with the intermediate structure of |32His97 at 3.7 ns (white), (c) Distance changes between the COM of 
the |3His97 sidechain and the COMs of aThe38 and aPro44 at the two switch regions plotted with x- 
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Figure 3 | % Distribution shift by high O2 partial pressure. Distributions 
of the intersubunit dihedral angle x calculated from 128 independent MD 
trajectories during 7-8 ns. The dihedral angles were sampled every 1 ps 
from 7 to 8 ns of each trajectory (in total 128,000 samples from 128 
trajectories). 

To obtain further structural insight, we investigated the x distri- 
butions of MD°2 and MD"^"' during 7-8 ns (Fig. 3). The two distri- 
bution curves were well fitted by the sum of three Gaussian functions 
as 

The fitted parameters (A/, dxv l^i) and the fitted curves are summar- 
ized in the Supplementary Table SI and Fig. S5. For MD'^''^ the main 
distribution gi{x) with the mean angle (Xi^) and width (5^1 ) of 
— 100.2 and 1.54 degrees was obtained. Meanwhile, for MD°^, the 
main distribution was —99.6 degrees mean angle and 1.75 degrees 
width. 

The comparison of MD°^ and MD'^''* distributions allows us to 
conclude that the solution environment of high O2 partial pressure 
enhances the T to R quaternary change from the following two 
aspects. First, from the 0.6 degrees peak position shift of the main 
distribution toward the R state (from -100.2 (MD^^O to -99.6 
(MD°^) degrees), the T state structure is biased toward the R state. 
Second, from the width increase from 1.54 (MD^'O to 1.75 (MD°^) 
degrees, the global structural fluctuation of the T state, which would 
contribute to the T to R quaternary change, is enhanced. The 
observed enhancement can be considered as the "non-site-specific" 
allosteric effect, or another allosteric regulation by O2 without bind- 
ing to the heme. 

Discussion 

Our ensemble MD simulation revealed that, even without the heme 
deformation by the O2 binding to the heme, the solution envir- 
onment of high O2 partial pressure itself enhances the quaternary 
structural change from T to R. This means that the traditional site- 
specific allosteric regulation by O2 binding to the heme is not neces- 
sarily the only one unique mechanism of O2 allostery and the 
additional non-site-specific regulation does exist. From a molecular 
point of view, we propose following two hypothetical mechanisms of 
the non-site-specific allostery of O2. 

First is the internal effect: O2 molecules enter the hydrophobic 
cavities inside HbA subunits^° and trigger the tertiary structural 
changes to enhance the T to R quaternary change. In MD°^ simula- 
tions, hydrophobic O2 molecules tended to escape from water and 
prefer the hydrophobic environment on the surface or inside HbA 
subunits^^. The O2 migrations between inside cavities would bring 
about the tertiary structural changes of HbA subunits since it was 




Figure 4 | Distribution of O2 oxygen atoms around switch region. 

Distribution of O2 oxygen atoms calculated from 128 MD trajectories 
during 7 to 8 ns in the vicinity of the switch region at the ai(32-interface 
switch region. Locations with density >0.05 A~^ and >0.1 A"^ are drawn 
with red wire frame and solid surface, respectively. Since the average 
density in the bulk solvent region (distance to the nearest protein atom 
>4.0 A) is about 0.005 A~^, the wire frame and solid surface regions are 10 
and 20 times denser than that in the bulk solvent, respectively. Four high 
density locations around the switch region are numbered in descending 
order of the peak density. 

observed that, in Mb, structural changes occur in response to the 
ligand migrations between cavities^ In particular, these reported 
structural changes in response to the CO transitions between the 
heme pocket and Xe4 cavity included displacements of the F helix, 
which forms the switch region of HbA P subunit^^ and hence would 
affect the stability of the T state structure. 

Second is the surface effect: the existence of O2 around the subunit 
interfaces weakens the contacts which stabilize the T state structure^ 
by affecting the surface residues directly and/or indirectly through 
perturbing the behavior of water molecules. In fact, there are 
obviously high O2 density locations in the vicinity of the switch 
region (Fig. 4). Thus it should be possible that the existance of O2 
directly affects the switch region residues through steric hinderance 
to weaken the contacts, while O2 also influences surface water mole- 
cules to indirectly affect the switch region residues. In general, the 
characteristic constants of water molecules near the protein surfaces, 
such as diffusion constant^^, can be different from that in bulk solv- 
ent. Moreover, it has been reported that the number of interfacial 
water molecules changes during quaternary changes in hemoglo- 
bins^^"^^, suggesting the important roles of the number of interfacial 
water molecules in the quaternary regulation. Thus, we can conjec- 
ture that the existence of O2 would perturb the behavior of water 
molecules, e.g., as the deviation of diffusion constants and the num- 
ber of interfacial water molecules around the surfaces, and, as a 
result, indirectly weaken the intersubunit contacts to enhance the 
T to R quaternary change. 

The O2 partial pressure applied in this work (0.55 mol/L) is about 
500 times higher than that in the O2 saturation concentration at 
ambient conditions (~1 mmol/L). We anticipate that the non-site- 
specific effects should be also observed under the ambient conditions 
for the following three reasons. First, our previous computational 
analysis with the same concentration 0.55 mol/L reproduced the rate 
constants of O2 entry into the binding sites of HbA subunits 
Assuming a kinetic model of O2 entry 

Hb^ + 02-^Hb^:02, (2) 

where x denotes subunit (x = a or p) and Hb^ : O2 is the situation of 
O2 in the binding site, the rate of O2 entry is A:^^^j.y[02] [Hb^] and 

the rate constants k^^^^ and A;fntry ^re calculated to be 45 and 99 
(|imol/L)"^s"\ respectively. These values are consistent with 
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experimentally observed values after temperature correction, 48-69 
and 81-131 (|imol/L)"^s"^ for a and p subunits, respectively. This 
consistency indicates that the concentration 0.55 mol/L is within the 
linearly extrapolatable range by the O2 concentration, as the rate 
constants can be estimated by the formula A;gj^^j.y[02][Hb^]. Second, 
the "effective" O2 concentration around HbA is higher than the bulk 
concentration. As discussed above, hydrophobic O2 molecules prefer 
the hydrophobic environment near the HbA surface and there are 
several high O2 density regions, whose density is an order of mag- 
nitude higher than that in the bulk solvent region (Fig. 4). Therefore, 
compared to the bulk concentration, the "effective" O2 concentra- 
tion around HbA is higher. This should partially narrow the concen- 
tration gap between the current simulation conditions and ambient 
conditions. Third, a preliminary ensemble MD simulation at one- 
tenth concentration, 0.055 mol/L (MD^^^^^ with 12 O2 and -12000 
H2O molecules) also revealed the non-site-specific effect. We exe- 
cuted 40 MD simulations for 8 ns at the 0.055 mol/L concentration 
and calculated the ^ distribution during 7-8 ns (Supplementary Fig. 
S6). The MD^°^°^ distribution is apparently shifted toward the R state 
as in MD^l 

From an experimental point of view, it is necessary to analyze the 
number of O2 molecules around or inside HbA subunits because the 
number of 02-bound heme in tetrameric HbA, which is the tra- 
ditional index of O2 allostery and is easily observed by spectroscopy, 
cannot capture the non-site-specific O2 allosteric contributions. In 
particular, as Tomita et al. discussed in Mb case^^, accurate measure- 
ment of the heme:02 stoichiometry with modern instrumentation is 
desired to verify the traditionally employed stoichiometry 1 : 1, which 
assumes that only one O2 binds to the heme and there is effectively no 
O2 in the hydrophobic cavities or near the surfaces of HbA subunits. 
The stoichiometry greater than 1 : 1 is plausible by hydrophobic 
interactions between O2 and HbA and means that HbA subunits 
can carry excess O2 in their cavities and/or surfaces, or the "effective" 
O2 concentration around HbA is higher than bulk concentarion, and 
support the existence of the non- site- specific effect. 

The current non-site-specific effect should provide a complement- 
ary mechanism to account for the cleavages of intersubunit contacts 
during the T to R quaternary change of HbA in high O2 concentra- 
tion environment from the molecular point of view. With respect to 
the displacements in the heme average strucutre, the 02-heme bind- 
ing itself brings about very small structural rearrangements: a com- 
parison of high resolution crystal structures of Mb revealed that the 
heme iron atom is 0.290 A displaced by ligand binding^^. Since the 
immediate displacements in the residues distant from the heme 
within 20 ps after ligand phosolysis are smaller than those neighbor- 
ing to the heme because of the elastic behavior of Mb^\ the imme- 
diate displacements propagated to the residues composing the 
intersubunit contacts must be smaller than the iron displacement. 
It is not obvious how such a small immediate displacements could 
cleavage the contacts, although it was in fact experimentally obserbed 
that the T to R quaternary change occurs after microseconds of 
ligand dissociation^^. Meanwhile, by the dynamic non-equilibrium 
MD approach in Mb, it was shown that a large structural fluctuation 
in the FG- corner (that of P subunit forms the switch region in HbA) 
was brought about by the ligand dissociation and recombination to 
the heme^\ Together with this structural fluctuation effect, O2 mole- 
cules located near the intersubunit contacts as in Fig. 4 can interact 
with the contacts, playing a complementary role in cleavage of the 
contacts. 

The non- site- specific effect can also be an important factor in 
allosteric structural regulation in other proteins. For example, in 
muscarinic receptors, orthosteric ligands were observed to function 
as an allosteric modulator^^ or weakly bind to the allosteric sites^^, 
indicating that the ligands can interact with multiple sites; not only 
the orthosteric site but also allosteric sites. In particular, for allosteric 
proteins in which the distance between orthosteric and allosteric sites 



is so far-away that direct structural perturbation between the sites 
seems to be unlikely, the multiple interaction sites between proteins 
and ligands, most of which are invisible by the X-ray crystallography 
because of their small occupancy, can help the structural rearrange- 
ments at the orthosteric sites. The concept of non- site- specific allos- 
tery should facilitate further understanding of allosteric regulation 
process depending on the concentration of effectors from the ato- 
mistic point of view. 

Methods 

The MD simulation force field parameters were identical to our previous work^^. All 
the MD simulations were executed by AMBER 9 pmemd module^"^. The initial HbA 
structure was retrieved from the T state deoxy-HbA crystal (PDB ID: 2DN2)^^ and 
immersed in the periodic boundary TIP3P water box with 12023 water molecules and 
two chloride counterions, yielding the MD"""* system. In addition, 120 water mole- 
cules were replaced by 120 O2 molecules to prepare MD*^^ system. We executed 128 
independent MD simulations for MD"''' and MD°^ systems, respectively, according to 
the following procedure with different initial atomic velocities. First, we executed 
300 ps high temperature NVT MD simulations at 700 K with the structure restraints 
on the protein. Second, equilibration NPT MD simulations for 200 ps at the target 
temperature 3 10 K and pressure 1 atm were executed with the protein restraint. Then 
production NPT MD simulations at 1 atm and 310 K were executed for 8 ns without 
any restraints. All the images of protein structures were drawn with VMD 1 .9.1^^. The 
density map in Fig. 4 was calculated by "VolMap Tool" in VMD with the parameters 
1.0 A resolution and 1.0 X radius atom size. 
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